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Abstract 

In this paper, with the corresponding formula for internal energy ob- 
tained in Ref. ^6. , combined with a simple calculation of the effective mass 
of interacting Bose particles, the behavior of the heat capacity of liquid 
*He is analyzed numerically for the entire temperature range. The results 
agree quite well with experimental data. 
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1 Introduction 

Notwithstanding a great number of papers (starting from Refs. [HE]) concerned 
with the microscopic study of Bose system's properties a good description of the 
heat capacity of hquid helium-4 in the whole temperature range has not yet been 
created. The first attempts were made by Brout [3] where it was shown in the 
first order of perturbation theory for the free energy of non-ideal Bose system 
that the presence of interaction does change the order of the phase transition. In 
[H [S] the thermodynamic functions of liquid helium-4 at all temperatures were 
obtained using the two-time temperature Green's function formalism. A good 
agreement of the specific heat at low temperatures with experimental data was 
obtained and the temperature of phase transition was calculated as Tc — 1.99 
K . The success of such an approach lies in the application of the experimentally 
measured structure factor of liquid helium-4 extrapolated to zero temperature [4] 
instead of the interparticle interaction potential. A good agreement of the heat 
capacity curve with relevant experimental data for the temperatures below the 
temperature of the A-transition was obtained in Ref. [B] where the calculations 
were made using the quantum-statistical approach based on the density matrix 
of Bose liquid. At a higher temperatures the specific heat curve was shifted 
upward almost in a parallel way. 

In Ref. [7] the density matrix formalism with the functional optimization 
of the Jastrow wave-function parameters was used to describe the properties of 
liquid helium. The results for internal energy agree well with experimental data 
for the temperatures below the critical one. It was shown that by taking into 
account the dynamic two-particle correlations only one can obtain the value of 



the critical temperature 3.4 K. Thus for good agreement with experiments one 
needs to take into account higher-order approximations, which are specificahy 
related to the concept of effective mass of the helium atom in a liquid. 

In recent years much attention has been paid to the study of the atom's 
effective mass in liquid helium because in this way part of the interaction could 
be taken into account accordingly to Feynman's idea [S]. However, there is no 
satisfactory formula for the effective mass of the helium atom in the liquid at 
arbitrary temperatures. Various scholars were mostly concerned with the value 
of the effective mass at T — > 0. Isihara and Samulski [S] have used the value of 
m* jm = 1.71 to agree the theoretically calculated sound branch of the excitation 
spectrum with the corresponding experimental data. In Ref. |10| the effective 
mass m*/m = 1.70 was obtained on the basis of the liquid helium-4 structure 
factor measurements. In Ref. [11] the interatomic potential was preserved as 
the input information, but in part the contribution of higher correlations was 
"transferred" to the kinetic energy term. In this way the mass of particles was 
renormalized that is somehow in correlation with the approach of Ref. [5] . As a 
result of such a renormalization the value of m* jm = 1.58 was obtained using 
Green's function method. It was shown in Refs. [l^ that the above-mentioned 
mass renormalization leads to the expressions obtained for the effective mass of 
the "^He impurity atom in liquid *He but with the replacement of the "pure" 
^He atom mass by the "^He atom mass. 

The aim of this paper is to calculate the heat capacity of liquid helium above 
the temperature of phase transition. The formula for the internal energy of Bose 
liquid, obtained in Ref. [6] with the help of the method proposed in Ref. [13] 
(where the effective mass is a free parameter of the theory), forms the basis 
of these calculations. Thus, further we calculate step by step with the help 
of thermodynamic perturbation theory the quasi-particle spectrum of the Bose 
system at the temperatures higher than the critical one, then we obtain the 
effective mass and numerically analyze the behavior of the heat capacity. 

2 Perturbation theory for the grand canonical 
potential at T > Tc 

Consider a collection of N spinless particles embedded into volume V . The 
Hamiltonian of the system which takes into account only pair interaction be- 
tween particles may be written using the secondary quantization language 

/Y2 1 AT 

k^^O kT^O 

Ho = ^{Sp - fJ.)a+ap. 
p 

The creation and destruction Cp operators of the particle with the momentum 
hp satisfy the usual bosonic commutation relations. The notations i'{k) stands 
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for the Fourier transform of the potential and Sp = h^p^ /2ra for the free-particle 
spectrum are introduced. It is more convenient to work in the grand canonical 
ensemble. That is why we introduced the chemical potential /x and fugacity 
z = e^^^{l3 = l/T, where T is the temperature) of the system that can be found 
with the help of the following equation 

Z(a+ap)=iV. (2) 
p 

Using the secondary quantization formalism it is easy to write down the Fourier 
transform of the particle density fluctuation operator 

Pk = -^Za+ap+k, (k^O). (3) 
p 

Further, our task is to calculate the partition function of a many-boson sys- 
tem above the temperature of phase transition. Of cause, the most interesting 
features of these calculations occur at the region in a close vicinity of the tem- 
perature of Bosc condensation. In the statistical operator, let us pass to the 
interaction representation and write down the partition function in the following 
way: 

{AT^ AT 1 

-l3—u{0) + l3—J2,.{k)\x 
k^O J 

X exp I - i J dr ^z/(fc)pk(T)p-k(r)|^ , (4) 



where the quantity 



Pk(T) = e^^Vkfi-^^" 



and the braces stand for statistical averaging with the Hamiltonian Hq. The 
first multiplier Zq is the partition function of the ideal Bose gas. The second one 
and the third one take into account the inter-particle interaction completely. 

Next, we rewrite, with the help of the Hubbard-Stratonovich transformation, 
the T-exponent in terms of the functional integral and we also average it by the 
states of the ideal Bose gas 

Z/Zo = exp z.(0)+/3— ZKfc)[ 

X J £)(^exp|-i^(l-Fz/(A;)n(a;„,fc))v?g</?-9 

(>3 qi qi J 

gi -I-.. .-1-91=0 
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where we use the notations q = (aj„, k) and aj„ = iirnT {n = 0; ±1; ±2; . . .) is 
the Matsubara frequency and 

E-EE- 

q tj„ 

Here the symbol J Dip denotes the integration over real and imaginary parts 
of the ipq variables from the half space of all possible values of q due to the 
symmetry if* — ip-q. The polarization operator 

HK, fc) ^ E G.{q')G,{q' + g) = ^ E ■ 7 ■ (6) 

Vp V ^ lUIn- Ek'+k + Ek' 

q' k' 

Here and thereafter is the filling factor of the ideal Bose gas. We also 
introduced notations for the symmetrical functions 

ni((7i, . ..,(?;) = E {Go{q)Go{q + qi)... Go{q + qi + . . . + qi-i) + permutations}. (7) 

The one-particle Green's function of noninteracting bosons is 

Goiq) - . 

lUJn - ek+ fJ. 

The representation Eq. ([5]) when the partition function is written in terms of 
functional integrals was used successfully in the theory of Fermi systems with 
Coulombic interaction in Ref. [14]. The fact that in the Gaussian approxi- 
mation of the calculation of the functional integral we recover Random Phase 
Approximation (RPA) correctly is a great advantage of our method. This ap- 
proximation in the case of interacting fermions generalizes at finite temperatures 
the well-known result of Gell-Mann-Brueckner for the high-density electron gas. 
The non-Gaussian part can be taken into account approximately by means of 
perturbation theory. 

The thermodynamic potential up to the first order of the perturbation theory 
(in two-sum approximation over the "4- vector" ) is 

= f^o + fii + f]2, (8) 

where the ideal gas contribution 

r!o = rEln(l-ze-^^''), (9) 
p 

and the RPA-part 

= — j.(o) + — E In |i + i.(A:)nK, fc)l - ^ E '^(^)- (10) 
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As we work in the grand canonical ensemble we take the average number of 
particles N of the function of the chemical potential /i. 

The one-loop contribution to the thermodynamic potential is 

^2 = 2 . Z^_VB2 n3(gi, ga, g3)n3(-gi, -92, -g3)/(gi)/(g2)/(g3) 

91 92 92 
91+92+93=0 

~ ^^7^XlXln4(gi,-gi,'72,-92)/('7i)/('?2)- (11) 

91 92 

The function f{q) = iy{k){ipqip^q) , and the correlator 

1 + i/(fc)n(a;„, k) 

were obtained in the Random Phase Approximation. The structure of expres- 
sions (flUl) . ((TT|) clearly shows that they are "in correlation" with the formulae 
obtained in Ref. [6] where an entirely different method of calculations was used. 



3 Renormalization of the one-particle spectrum 



It is clear that the basis of further calculations is fully determined by the 
renormalization of the quasi-particle spectrum. For our analysis we use RPA. 
Notwithstanding the simplicity of this approximation it "catches" certain im- 
portant features of the behavior of the system. It is not surprising because RPA 
effectively sums up an infinite set of terms of the perturbation theory divergent 
near phase transition point. 

At first let us use the thermodynamic equality —dil/dfi — N to find the 
average number of particles in the system. The explicit calculation of the cor- 
responding derivative with the first two terms of Eq. ([5]) gives 



dup 



1 u[k) 



2V(3 



l + u{k)W{uJn, 



k) 



Let us construct the Bose filling factor with a new spectrum using the expression 
in braces. Making use of equality dup/dpL = —dup/dsp we finally obtained the 
formula for the renormalized one-particle spectrum 



where the correction to the quasi-particle spectrum is 

Ae 1 '^(fc) [ 1 

pV ^ 1 + iy{k)U{ujn,k) { e|k-p| - Ep - iu 



Sk - lljJn 



(13) 



(14) 
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and the correction to the chemical potential is 



N 
V 



1 



Kfc) 



i/(fc)n(cj„,fc) el 



1 



MO 



(15) 



It is easy to obtain the above-mentioned expression for the spectrum in a differ- 
ent way. To do this one has to recall that variational derivative of the fi-potential 
with respect to £p equals the renormalized one-particle filling factor n*. After 
simple calculations we obtained the following formula: 



^.(0)-^"^ 



V 



de-r, 



1 

W 



2VI3 ^ 1 

9 



£|k-p| 



j/(fc)n(w„, k) 

V (w„ 



+ 



"-Ik-pl 



(e|k-p| 



+ (w„ 



-w„) 

-l^ri) 



drip 



(16) 



It is easy to argue by making summation over the wave- vector p of the left-hand 
and right-hand sides of the previous equality that the second term in braces 
vanishes. So, after the summation of expression (jl6|) we arrive at equality (jl2p 
and thus we get the same expression for the correction to a one-particle 
spectrum again. It is interesting to note that the calculation of the variation 
derivative Sfl/Sup in RPA gives the same result. Finally, Eq. (fT4)) coincides 
with the result derived in Ref. jl5| where calculations were made in terms of 
temperature Green's function technique. 

Let us analyze expression (fH|) . First, the potential problems with the inte- 
gration over the wave-vector may occur only in the critical region and at the 
zero frequency a;„. That is why we write down this term apart and immediately 
set apart the Hartree-Fock-like term 



Ae, 



1 



£|k-p| 



- £r, 



1 

£k 



1 ^ v{k) 
Wi^,l + HkMk) 

^^^^ I '^(/^^Ik-pl - Psp) - n{P£k) 

1 ^ t/2(fc)n(c^„,fc) 
pv ^^i + v{k)n{uj^,i 



|k-p| 



- £,1 



1 



Pek] 



Mk-p| 



Ek - liOn 



(17) 



where H{k) = 11(0, fc), n{x) ~ l/(e^ — 1) and p = N/V is equilibrium density 
of the system. Secondly, to go further we have to investigate the properties of 
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polarization operator ^ 



n(w„, k) 



(27r)2 2k 



dx ■ 



In 



(k/ko)'^ + 2xk/ko — iur, 



(fc/fco)^ — 2.xfc/fco — iUr, 



(18) 



here for convenience the fohowing notations are used: fco — \/2mT/h, u„ = 27rn. 
We are interested in a long wave-length behavior of the polarization operator. 
To find the leading order asymptote of n(w„, k) at zero frequency in the critical 
point it is sufRcient to replace the Bose filing factor l/(e^ — 1) in integral (dHJ 
by Then after a simple integral calculation we get n(fc — > 0) = jSk^K^k). 



For higher temperatures (/i ^ 0) 



g-mief'n) /53/2(e^^) + 



(19) 



The dots stand for higher than quadratic terms in the expansion over the 
wave- vector. It is easy to see from definition © that for non-zero frequencies 

n(w„,o) = o. 

Now, let us consider the contribution to the one-particle spectrum (|14p from 
zero frequency. Using designation for this term of spectrum it is easy to 
write down 



Ae;^ = (fco/7r)2p / dx 



v{2xp) 



1 + pv{2xp)Tl{2xp) 



"^In 


x + 


1 








- 1 


2 


X — 


1 





(20) 



For self-consistency of our calculations, especially near the critical point, the 
chemical potential fi should be changed by p* in the right-hand side of Eq. 
(PD|) (the critical point is determined by the equation yu* = 0, respectively). 
Admittedly, the ideal gas dispersion relation should be replaced by the exact one- 
particle spectrum, but further analysis will not be influenced by this replacement 
qualitatively. 

Let us consider the value of the integral in Eq. (f20| at a small p and assume 
for definiteness that the temperature is higher than the critical one. Then 
substituting Ii{2xp) — > 11(0), and v{2xp) i/(G) we obtain 



As; = {ko/TrYp- 



p^(o)n(o) 



j dx 


■^In 


x + 


1 










- 1 




2 


X — 


1 





It turns out that this integral equals zero identically. Moreover even after the 
substitution of n(fc) — > (|19p . the integral in Eq. (pn)) equals zero too. Thus it 
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is shown, with a reahstic restriction on the Fourier transform of the potential 
energy, i.e. the absence of hnear and quadratic terms in the expansion of I'ik) 
at a small fc, that Ae^ — o{p'^). 

The situation is quite different in the critical region. Here the leading order 
asymptote is Ae^ ~ ln(p) (it is not hard to ascertain taking into account the 
properties of n(fc) ^ 1/fc in this region), which obviously is a hint at the follow- 
ing behavior of the one-particle spectrum Ae^ ~ p^^"^ {rj — 0) at the critical 
temperature. Clearly one cannot obtain this result using a simple perturbative 
approach. 

Hence, having separated the non-analytic problematic part of the spectrum 
(the second term in Eq. p7)) ) we can consider the "non- universal" one, i.e. the 
remainder of Eq. (ITTI) . Precisely this expression will determine the observed 
non- universal properties of the Bose liquid. The latter calculations are linked 
to the summation over the Matsubara frequency in the last term of Eq. (|17p 
and coincide with those in Ref . [15] . That is why we do not dwell on the details 
of these calculations. Now the leading-order non- vanishing term of the quasi- 
particle dispersion relation is quadratic over the wave-vector. We recall that 
ASp = o(p^) and hence its contribution is not significant. So, for reasons of 
simplification we assume the spectrum to be a quadratic free-particle one, but 
with the normalized mass. 

As we single out the "problematic" contribution to the quasi-particle spec- 
trum expanding second and third sums of Eq. (jl7l) into a series in p we obtain 
for the one-particle spectrum 

^; = As," + ^. (21) 



The effective mass is 



where the quantity 



m/m* = 1 - A(T), (22) 



= 4E """" 



3N 



^|4±|[n(/3£fc)-l//3e.] 



ak{al - 1) 

+ 2[l//3efc-/3£fen(/3efc)[l + n(/3£fc)]] }• (23) 

Here we use the following notations: E/- = ak£k is Bogoliubov spectrum, and 
ak = a/1 -|- 2pv[k)/ek- The effective mass in the low-temperature region is 
always larger than its "bare" one which means that the renormalized tempera- 
ture of the Bose condensation of interacting particles is always lower than the 
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critical temperature of the ideal gas. This is the most important result of Eqs. 
([22]) and ([23]). At low temperatures carefully calculating the limit of /3 — >■ it 
is easy to ascertain that the effective mass tends to the mass of particles. It 
is important that the temperature-independent part of formula (j23p coincides 
with the effective mass derived in Ref. |12| where a different method was used 
for the calculations. 

At the end of this section one more remark on the applicability of the formula 
for A(T) has to be made. The approximation of the exact spectrum of collective 
modes by the Bogoliubov spectrum on the one hand made it possible to obtain 
the analytical expression (|^^ . on the other hand it brought us beyond the 
limits of RPA. It is hard to assess the accuracy of such a trick, but it can be 
justified only considering that the Fourier transform of a two-particle potential 
is a rapidly decreasing function. Then the main contribution to the integral over 
the wave-vector comes from a lower limit of integration where the Bogoliubov 
spectrum completely coincides with the exact one. 



4 Internal energy and heat capacity 

The expression for the grand canonical potential derived in the second section 
of this paper is applicable only for the temperatures higher than the critical 
one. To describe the A-transition phenomenon, in particular thermodynamic 
functions, let us use an approach based on the density matrix of the Bose liquid 
[13] . Using this approach the dependence of an internal energy of the Bose liquid 
on the effective mass of the helium-4 atom was found in the approximation of 
pair-particle correlations in Ref. [5]. In the case of m* /m = 1 in Ref. [B] a good 
agreement of the heat capacity with experimental data in the region below the 
critical point was obtained, but above the temperature of phase transition the 
heat capacity curve was shifted upward. Hence, the calculation of heat capacity 
for the case of m*/m ^ 1 is an interesting problem. 

We take the expression for the internal energy in the approximation of pair- 
particle correlations from Ref. 
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2 7,2 



1 



2m 



4 ^ 

-E 



+ al - 1) S{k) + 



2l2 



H'k 



16 ^ 2m 

k#0 



1 m* 

-E 

2 ^ 2m 



k^^O 



1 



dSoik) 



Afc5o(fc) 



dp 



sinh {f3Ek) 



sinh (/3e* 

(24) 



Here the following notations are introduced: 



^o(fc) = l + 2^n^ + l 



^^P^^|p+k|> 



(25) 



p+k^^O 
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is the structure factor of the ideal Bose gas with a renormaUzed mass. 

The quantity no is an average number of particles of the ideal Bose gas with 
zero momentum and zq is its fugacity. 

The next quantity 

is the pair structure factor of the Bose liquid, and 

Afc = ak tanh {/3Ek/2) - tanh il3el/2) . (27) 

Obviously, if we turn off interparticle interaction at = 1, m* — m and take 
into account that sound velocity in the ideal Bose gas at T = equals zero c = 
(there are no zero density fluctuations in the Bose system) and Afc = (|27|) we 
obtain a well-known formula for the energy of the ideal Bose gas 

Let us analyze the total energy Eq. p4|) in the low temperature region where 
it coincides with the formulas obtained in Ref. [13l [16] . At low temperatures 
T — >■ 0, when only small values of the wave-vector k are important in the 
expression for the spectrum Ek wc obtain: 



where the ground state energy is 

mc^ ly-^h^P . ,2 1 ^K^k^ 1 ( 1 

^0 = N— {ak^lf + -y,-^ "fc — 

2 4 ^-^ 2m 16 2m \ afc 

and the heat capacity, respectively. 



Exactly the same temperature dependence of the heat capacity of liquid 
'^He at T — > is observed. Since we obtained the correct behavior of the heat 



capacity at low temperatures using the energy of the Bose liquid Eq. (|24p 
we expect to derive the correct behavior of the heat capacity for the entire 
temperature range. 

5 Numerical results 

Our numerical calculations are carried out at the equilibrium density of liquid 
helium p = 0.02185 A^^^ mass of particles m = 4.0026 u, sound velocity c = 
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238.2 m/s in the hmit of T -> [T^, and at the critical temperature of the ideal 
gas Tc = 3.138 K. We use the liquid structure factor extrapolated to T = [4] 
as the output information, instead of the interparticle potential, i.e. 

- -S^y (29) 

where S'^^P{k) is the experimentally measured structure factor at T = 0. 

It is logical to start calculations with the formula for the renormalized one- 
particle spectrum and thus with the formula for the effective mass of particles. 
Despite the complexity of the last sum over the wave- vector in Eq. the main 
contribution to the effective mass arises from the temperature- independent part. 
In Fig. 1 the dependence of a dimensionless value of the effective mass m* /m as 
a function of temperature is presented. Formally we extrapolated a curve of the 
effective mass in the condensate region where obviously it becomes a parameter 
of the theory. It is important that at zero temperature m* coincides with the 
effective mass of the impurity atom in the Bose liquid [T2] . 

The first three terms of a low-temperature expansion are (T ^ mc^, T <^ 

A(T ^ 0) = Ao + AiT - Ag/aT^/^ -I- o{T^''^). (30) 
where A-coefficients are 

Ao^^v ("r^)^,^o.4i, 

^^-4E4^-0.31, 



^^/^ - 3^(^/2) (2^) /^ = «-2 
Then the effective mass equals approximately 

^0.59-0.31T + 0.24r3/2 



This formula reproduces the curve in Fig. 1 quite well up to the critical tem- 
perature. 

Now we are in position to calculate the renormalized temperature of the Bose 
condensation. We can find Tc using condition n* = iV at the zero value of 
the renormalized chemical potential /i*. A simple calculation gives Tc = 1.94 K 
that agrees quite well with experimental measurements of the temperature of 
the A-transition Tc^p = 2.17 K despite the simplicity of the approximations. 

Let us pass on to the heat capacity calculation: 
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m*/m 

2.0- 




Figure 1: Temperature dependence of the fraction m* /m 

We calculate heat capacity using the difference method and build the plot 
of its temperature dependence T/Tc- 

A comparison of different heat capacity curves is depicted in Fig. 2. As 
is seen from the comparison of the calculated curve 1 with the experimental 
one the agreement is quite good at low temperatures (0 < T/Tc < 1). At the 
temperatures T/Tc> 1 the inconsistency occurs: the behavior of the calculated 
heat capacity is very similar to the behavior of the experimental curve, but 
shifted upward almost in a parallel way. This inconsistency is related to the 
fact that three- and four-particle correlations should be taken into account for 
the quantitative description. The contribution of three- and four-particle cor- 
relations, as is shown in Refs. [16j[T8], improves significantly the ground-state 
results and gives a fairly good agreement at T — 0. 

Further, let us calculate the heat capacity of liquid helium-4 with taking 
into account the effective mass of the helium atom in the liquid (curve 2). At 
low temperatures the heat capacity with taking into account the effective mass 
of the Bose particles practically coincides with curve 1, which agrees well with 
experimental data. This shows a weak dependence of the heat capacity on the 
effective mass below the temperature of phase transition. As is seen from Fig. 
2, the calculated curve 2 (unlike curve 1) agrees quite well with the experimen- 
tal one. It is related to the fact that by using the effective mass we partially 
take into account a contribution from three- and four-particle correlations. It is 
not surprising that in close vicinity of the Bose condensation point the theoreti- 
cally calculated heat capacity deviates most significantly from the experimental 
curve. It is solely related to the inconsistency of our description near the critical 
point because the non-analytical part (PP)) of the one-particle spectrum, which 
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Figure 2: Heat capacity of liquid heliuni-4. Curve 1 is the calculated heat capac- 
ity in the pair-correlation approximation (j24p [6]; curve 2 is the heat capacity 
in pair-correlation approximation (j24p taking into account the effective mass of 
the helium atom in liquid ()22|) . The circles show experimental data from Refs. 
[19, 20, 214. 



makes a significant contribution in the thermodynamic functions at T — > Tj,, is 
disregarded in our approach. One has to use renormalization group methods 
[231 HH di] for the correct description of the heat capacity in this temperature 
region. 



6 Conclusions 

In this paper we succeed in deriving quite well an agreement of the heat capacity 
curve of liquid helium with experimental data practically for all temperatures. 
The application of the thermodynamic functions of the Bose liquid obtained 
with the help of the Hamiltonian averaging combined with a one-particle spec- 
trum derivation were the key moments of the calculations. Notwithstanding the 
simplicity of the spectrum calculation we obtained quite interesting results. In 
particular, we can decompose the part of the quasiparticle spectrum that is re- 
sponsible for the non-analyticity in the Bose condensation point and show that 
this term of spectrum has no effect on the physical observables in the undercrit- 
ical temperature region. So, an attempt is made to justify microscopically the 
idea that the A-transition in a real quantum liquid is very similar to the Bose- 
Einstein condensation phenomena of the ideal gas "slightly" deformed by the 
interaction between the particles (keeping in mind the non-universal properties 
of the system). 
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The calculation found that the long wave-length limit of the "non-universal" 
part of the one-particle spectrum is quadratic over the wave-vector, i.e. very 
similar to the dispersion relation of the ideal gas but with a new mass. In the 
general case it is shown that this new mass at low temperatures is always greater 
than the mass of particles, and thus, the presence of the interaction at least in 
our approximation always lowers the critical temperature. 

Another feature of the developed theory, perhaps a bit unexpected, is that 
even this simple temperature dependence of the effective mass improves the 
behavior of the heat capacity curve in the undercritical region and does not affect 
it in the condensate phase. Hence, the quantum-statistical approach based on 
the density matrix is suitable for describing thermodynamic properties of such 
a strongly-interacting Bose liquid as the helium-4 liquid not only in the limits 
of low and high temperatures, but for the entire temperature range. 
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